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Abstract 

We give an 0(n 15 logn) algorithm that, given a directed planar graph 
i— i with arc capacities, a set of source nodes and a single sink node, finds a 

m maximum flow from the sources to the sink . This is the first subquadratic- 

time strongly polynomial algorithm for the problem. 

lHj 1 Introduction 

The study of maximum flow in planar graphs has a long history. In 1956, Ford 
and Fulkerson introduced the max si-flow problem, gave a generic augmenting- 
path algorithm, and also gave a particular augmenting-path algorithm for the 
£f) case of a planar graph where s and t were on the same face (that face is tra- 

ditionally designated to be the infinite face). Researchers have since published 
qq many algorithmic results proving running-time bounds on max si-flow for (a) 

(~^) planar graphs where s and t are on the same face, (b) undirected planar graphs 

where s and t are arbitrary, and (c) directed planar graphs where s and t are 
arbitrary. The best bounds known are (a) 0(n) [IT], (b) O(nlogn) [BJ, and (c) 
O(nlogn) [I]. 

Schrijver [JJ5] has written about the history of this problem. Ford and Fulk- 
erson, who worked at RAND, were apparently motivated by classified work of 
Harris and Ross on interdiction of the Soviet railroad system. (Of course, Har- 
ris and Ross were interested in the min cut, not the max flow, as seems to 
be true for most applications.) This article was downgraded to unclassified in 
1999. It contains a diagram of a network that models the Soviet railroad system 
indicating "ORIGINS" (sources) and what is apparently a sink (marked "EG"). 

In max-flow applied to general graphs, multiple sources presents no problem: 
one can reduce the problem to the single-source case by introducing an artificial 
source and connecting it to all the sources. However, as Miller and Naor [T?] 
pointed out, this reduction violates planarity unless all the sources are on the 
same face to begin with. Miller and Naor raise the question of computing 
a maximum flow in a planar graph with multiple sources and multiple sinks. 
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Figure 1: The soviet rail network 
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Even when there is one sink, until now the best known algorithm for computing 
multiple-source max-flow in a planar graph is to use the reduction in conjunction 
with a max-flow algorithm for general graphs. That is, no planarity-exploiting 
algorithm was known for the problem. 

There are workarounds. For example, in the Soviet rail network, there are 
two faces that together include all the sources, so solving the instance can be 
reduced to solving two single-source max flows in a planar graph. However, a 
more realistic motivation comes from selecting multiple nonoverlapping regions 
in a planar structure. 

Consider, for example, the following image-segmentation problem. A grid is 
given in which each vertex represents a pixel, and edges connect orthogonally 
adjacent pixels. Each edge is assigned a cost such that the edge between two 
similar pixels has higher cost than that between two very different pixels. In 
addition, each pixel is assigned a weight. High weight reflects a high likelihood 
that the pixel belongs to the foreground; a low-weight pixel is more likely to 
belong to the background. 

The goal is to find a partition of the pixels into foreground and background 
to minimize the sum 

weight of background pixels 

+ cost of edges between foreground pixels and background pixels 

subject to the constraints that, for each component K of foreground pixels, the 
boundary of K forms a closed curve in the planar dual that surrounds all of K 
(essentially that the component is simply connected). 

This problem can be reduced to multiple-source, single-sink max-flow in a 
planar graph (in fact, essentially the grid). For each pixel vertex v, a new vertex 
v' , designated a source, is introduced and connected only to v. Then the sink 
is connected to the pixels at the outer boundary of the grid. 

New result We prove the following: 

Theorem 1.1 There is an 0(n l ' b \ogn) algorithm to compute multiple- source, 
single-sink max flow in an n-node directed planar graph. 

Before our work, the best strongly polynomial bound for the problem is 
0(n 2 logn), which comes from the reduction to general graphs and then use of 
an algorithm such as that of Goldberg and Tarjan [§]. For integer capacities 
less than U, one could instead use the algorithm of Goldberg and Rao 0, which 
leads to a running time of 0(n 15 logn log U), or the planarity-exploiting min- 
cost flow algorithm of [T2], which gives a bound of 0(n 1595 log U) that depends 
on fast matrix multiplication and interior-point methods. However, even if one 
assumes integer capacities and U — Q(n), our planarity-exploiting algorithm is 
asymptotically faster. 

We have learned (personal communication) that Borradaile and Wulff-Nilsen 
have independently proved the same theorem. 
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1.1 Organization 

The structure of the paper is as follows. In Section [2] we give some definitions 
and general technical background. In Section [3j we give the main algorithm. 
Finally, in Section [4] we describe how to efficiently convert a feasible preflow to 
a feasible flow. 

2 Preliminaries 

2.1 Embedded Planar Graphs 

A planar embedding of a graph assigns each node to a distinct point onto the 
sphere, and assigns each edge to a simple arc between the points corresponding 
to its endpoints, with the property that no arc-arc or arc-point intersections 
occur except for those corresponding to edge-node incidence in the graph. A 
graph is planar if it has a planar embedding. 

Assume the graph is connected, and consider the set of points on the sphere 
that are not assigned to any node or edge; each connected component of this 
set is a face of the embedding. 

It is convenient to designate one face as the infinite face (by analogy to 
embeddings on the plane). With respect to a choice of the infinite face, we say 
a Jordan curve strictly encloses an edge or node if the Jordan curve separates 
the edge or node from the infinite face. Similarly, for a subgraph, the choice 
of infinite face foo for the whole graph induces a choice of infinite face for 
each connected component of the subgraph, namely that face of the connected 
component that contains /oo. 

In implementations, an embedding onto the sphere can be represented com- 
binatorially, using a rotation system. 

2.2 Flow 

Let G be a directed graph with arc set A, node set V and sink t. For notational 
simplicity, we assume here and henceforth that G has no parallel edges and no 
self-loops. 

We associate with each arc a two darts d and d', one in the direction of a and 
the other in the opposite direction. We say that those two darts are reverses of 
each other, and write d = rev(d'). 

A flow assignment /(•) is a real- valued function on darts that satisfies anti- 
symmetry: 

/(rev(d)) = -/(d) (1) 

A capacity assignment c(-) is a function from darts to real numbers. A flow 
assignment /(•) respects capacities if, for every dart d, f(d) < c{d). Note that, 
by antisymmetry, /(d) < c(d) implies /(rev(d)) > — c(rev(d)). Thus a negative 
capacity on a dart acts as a lower bound on the flow on the reverse dart. In this 
paper, we assume all capacities are nonnegative, and therefore the all-zeroes 
flow respects the capacities. 
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For a given flow assignment /(•), the net inflow (or just inflow) node v 

is inflow/ («) = EoeAiead(a)=« /(°) ~ E ae A:tail(a)=«- The out fl ow of u is 
outflow/^) = — inflow/ (w). The waZwe of /(•) is the inflow at the sink, inflow/(i). 

A flow assignment /(•) is said to obey conservation if for every node v other 
than f, outflow/ (v) > 0. 

A supply assignment a(-) is a function from the non-sink nodes to MU {oo}. 
For any node v, a(v) specifies the amount of flow that can originate at v. A 
flow assignment /(•) is said to respect the supplies cr(-) if, for every node v other 
than the sink t, outflow/(t;) < cr(v). In this paper, we assume all supply values 
are nonnegative. 

A flow assignment is a feasible preflow if it respects both capacities and 
supplies. A feasible preflow is called a feasible flow if in addition it obeys 
conservation. In this paper, we give an algorithm to find a maximum (feasible) 
preflow, and then an algorithm to convert that preflow to a maximum (feasible) 
flow. 

The residual graph of G with respect to a flow assignment /(•) is the graph 
G f with the same arc-set, node-set and sink, and with capacity assignment c/(-) 
and supply assignment c/(-) defined as follows: 

• For every dart d, Cf(d) = c(d) — f(d). 

• For every node v, o~f(v) — o~{v) — outflow/(i;). 

Single-source limited max flow 

For a particular node s, a limited max st-flow is a flow assignment /(•) of 
maximum value that obeys capacities and for which inflow/ (v) = for every 
node except s and t and such that outflow/ (s) < o~(s). An algorithm for ordinary 
max st-flow can be used to compute limited max si-flow by introducing an 
artificial node s' and an arc s's of capacity c(s), and running the algorithm on 
the transformed graph. This transformation preserves planarity. Since there 
is an 0{n log n) algorithm for max si-flow in a planar directed graph pQ, we 
assume a subroutine for limited max si-flow. 

2.3 Jordan Separators for Embedded Planar Graphs 

For an n-node planar embedded simple graph G, we define a Jordan separator 
to be a Jordan curve S such that, for any arc a of G, the set of points in the 
sphere corresponding to a either (i) does not intersect S or (ii) coincides with 
a subcurve of S. We require in addition that, if the two endpoints of a are 
consecutive nodes on S, then (ii) must hold. The boundary nodes of S are the 
nodes S goes through. 

We say a Jordan separator is balanced if at most 2n/3 nodes are strictly 
enclosed by the curve and at most 2n/3 nodes are not enclosed. 

Miller |16j gave a linear-time algorithm that, given a triangulated two- 
connected n-node planar embedded graph, finds a simple cycle in the graph, 
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consisting of at most 2y/2^/n nodes, such that at most 2n/3 nodes are strictly 
enclosed by the cycle, and at most 2n/3 nodes are not enclosed. 

To find a balanced Jordan separator in a graph that is not necessarily tri- 
angulated or two-connected, add artificial edges to triangulate the graph and 
make it two-connected. Now apply Miller's algorithm to find a simple cycle sep- 
arator with the desired property. Viewed as a curve in the sphere, the resulting 
separator S satisfies the requirements of a balanced Jordan separator, and it 
has at most 1\/2yjn boundary nodes. 



3 The Algorithm 

The main algorithm finds a maximum preflow in the following, slightly more 
general, setting. 

• Input: 

— A directed planar embedded graph G, 

— a sink node t, 

— a nonnegative capacity assignment c(-), and 

— a nonnegative supply assignment &(■). 

• Output: A feasible preflow G of maximum value. 

We present the main algorithm as a recursive procedure with calls to a single- 
source limited-max-flow subroutine. We omit discussion of the base case of the 
recursion (the case where the graph size is smaller than a certain constant.) Each 
of the recursive calls operates on a subgraph of the original input graph. We 
assume one global flow assignment /(•) for the original input graph, one global 
capacity assignment c(-), and one global supply assignment cr(-). Whenever the 
single-source limited-max-flow subroutine is called, it takes as part of its input 

• the current residual capacity function c/(-) and 

• the current residual supply function cr/(-)- 

It computes a limited max flow assignment /(•), and then updates the global 
flow assignment /(•) by f{d) :— /(d) + f(d) for every dart in the subgraph. 
In the pseudocode, we do not explicitly mention /(•), c(-) cr(-), c/(-), or &/(•)■ 
The pseudocode for the algorithm is given below. We assume that the sink 
is on the boundary of the face designated the infinite face. 

The algorithm proceeds in iterations as long as the current graph, G,;, con- 
sists of more than Nq nodes, for some constant Nq to be specified later. For 
graphs of constant size, output the solution in constant time. At iteration i it 



finds a small Jordan separator Si in Gj_i as described in Section 2.3 Let Hi 
be the subgraph of G,_i enclosed by Si. Intuitively, one would like to think of 
Si as the external face of Hi. However, Si might cross some earlier Sj, so Si 
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Algorithm 1 MultipleSourceMaxPreFlow(graph Go, sink t) 



I: triangulate Go with zero-capacity edges. 

2: i := 

3: while Gi consists of more than N nodes do 

4: i := i + 1 

5: find a Jordan separator Si in Gi_i 

6: let Hi be the subgraph of Gj_i enclosed by Si 

7: let Gi be the external face of Hi. 

8: let Bi be the set of cycles {Cj : Cj is contained in Hi} 

9: for G in i?i do 

10: designate one of the nodes of G as an artificial sink t' and add artificial 

infinite-capacity edges parallel to G 

li: MultipleSourceMaxPreFlow(i/i, f ) 

12: remove the infinite-capacity artificial edges 

13: let Gi be the subgraph of Gj_i that is not strictly enclosed by Si 

14: for G in {Cj} do 

15: for every node v of G do 

16: limited max-flow from i; to t in G 



does not entirely coincide with edges of Hi (recall that the separator procedure 
first triangulates the input). To overcome this technical issue we consider Gj, 
the external face of Hi, instead of just Si. 

For every cycle Cj that is contained in Hi we compute a maximum prcflow 
to Cj in Hi. We then set Gi to be the part of Gj_i not strictly enclosed by Si 
and continue to the next iteration. 

When all iterations are done, for every node v of every cycle Cj we compute 
a maximum v-to-t flow in Gq. However, we do not push more than o~(v), the 



excess flow present at v. We call this step limited max-flow (Line 16) 



3.1 Correctness of Algorithm [T] 

Definition 3.1 (Admissible path) A u-to-v path P is called admissible ifo~(u) > 
and if P is residual. 

Lemma 3.2 Fix an iteration i of the while loop. At any time in that iteration 



after Line 11 is executed for some cycle G', there are no admissible to-C paths 
is Hi . 

Proof: By induction on the number of iterations of the loop in Line H For 



the base case, immediately after Line 11 is executed for cycle G', by maximality 



of the preflow pushed when the edges of C had infinite capacity, the lemma 
holds. Assume the lemma holds before Line 11 is executed for cycle C" and 



let / be the flow pushed in that execution. Assume for contradiction that after 
the execution there exists an admissible u-to-C' path P in Hi for some node 
u e Hi. 
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Figure 2: A possible situation in the proof of Lemma [3. 2| P is shown in solid 
blue, Q in dashed red. 



If P was residual before / is pushed, then a(u) must have been zero at that 
time. Since P is admissible after he push, a(u) > after / is pushed. Therefore, 
before the execution, there must have been an admissible x-to-u path R in Hi 
for some x £ Hi. Thus, R o P is an admissible x-to-C path in Hi before the 
execution, a contradiction. 

If P was not residual before / was pushed, there must be some dart of P 
whose reverse is used by /. Let d be the latest such dart in P. The fact that 
rev(d) is assigned positive flow by / implies that before / is pushed there exists 
an admissible path Q from some node x € Hi to head(d), see Fig. [2} By choice 
of d this implies that QoP[hea.d(d),v] is an admissible x-to-C path in Hi before 



Line 11 is executed for cycle C", a contradiction. QED 



Lemma 3.3 Just before the loop in Line \14\ is executed, for every i, there are 
no admissible to-Bi paths in Hi. 



Proof: Lemma [3.2| implies that, for every i, at the end of iteration i of the 
while loop, there are no admissible to Bi paths in Hi. 

Since for j > i Hi n Hj C d and since Ci S Bi, there are no admissible 
to-Bi paths in Hi at any later iteration as well. QED 



Lemma 3.4 Just before the loop in Line 14 is executed for the first time, there 
are no v-to-t admissible paths for any node v G Go — IJ^- Cj . 

Proof: Let v be a node of Go that does not belong to any Cj. Let i be the 
unique index such that v £ Hi. Observe that any v-to-t admissible flow path 
in Gq must visit some node of Bi before getting to t, so it consists of a v-to-C 



admissible path in Hi for some C € Bi, contradicting Lemma 3.3 QED 



Lemma 3.5 For any node u, if there are no u-to-t admissible paths before an 
execution of Line then there are none after the execution as well. 
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Proof: If a(u) = before the execution then o~(u) = after the execution as 
well, so there are no admissible u-to-t paths. 

Otherwise, there is no u-to-T residual path before the execution, let / be 
the v-to-t flow pushed at Line[l6j Assume for contradiction that after the ex- 
ecution there exists a u-to-t residual path P. Since P was not residual before 
/ was pushed there must be some dart of P whose reverse is used by /. Let 
d be the earliest such dart in P. The fact that rev(d) is assigned positive flow 
by / implies that before / is pushed there exists a residual path Q from tail(d) 
to t. By choice of d this implies that P[tail(d)] o Q is a residual (and therefore 
admissible) u-to-t path before the execution, a contradiction. QED 

We can now prove the correctness of the algorithm claimed in Theorem 



mimc- 



Proof: (Of correctness of algorithm in Theorem |1.1[ ) By Lemma 3.4 
diately before Linc |14| there are no u-to-t admissible paths in Go for any node 



v G G 
terminates 
for v G (J 



By Lemma 



3.5 



there are no such paths after the loop in Line 



14 



Since the executions of Line [161 eliminate all v-to-t admissible paths 



• Cj , there are no admissible paths to t in Go upon termination, 
the flow computed is a maximum preflow. QED 



so 



3.2 Running Time of Algorithm [T] 

Lemma 3.6 Every cycle Cj appears at most twice as the cycle stored by the 
variable C in the loop in Zme[9| 

Proof: Consider Cj. It appears as the cycle stored by the variable G in the 
following two cases: 

1. when i = j (i.e., when Cj is the external face of Hi). 

2. when Cj is contained by some Hi but is not the external face of Hi. 

Note that case (2) implies that Cj has some dart that is strictly enclosed by Si, 
so this can happen for exactly one value of i since the subgraph strictly enclosed 
by Si is not part of G^+i. Thus, Cj is not contained by any H- with i' > i. QED 

We first consider the cost of the recursive calls in Line [TT1 In both cases in 



Lemma 3.6 the recursive call is on the graph Hi, so if T(n) denotes the running 
time of Algorithm^ on an input graph with n nodes, the cycle Cj contributes 
at most T(\Hj\) + T(\H p ^\), where p(j) is the unique value such that Cj is 
contained by H p ^ but is not the external face of H p (j\. Therefore, the total 
time required by all recursive calls is 

Y,T(\Hj\)+T(\H p(j) \). 
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Observe that, since for every i, < \Hi\ — 2\/2|Gj_i| < ||Gj_i|, we 

have \Hi\ < \Hj\ for i > j. Also note that, if Cj is not the external face of Hi 
and Hi contains Cj, then i > j. Therefore, ^2jT(\Hj\) + T(\H p ^\) is bounded 
^ E 3 T(\H,\) +T(\H J+1 \) =T(\H 1 \) +2J2^ 2 T(\H 3 \). 

Lemma 3.7 Y,% kl l^f' 5 < ^\G kl -i\^ ■ 

Proof: By induction on ki — k\. Recall that, for every j, \Hj\ < ||Gj_i| + 
2V2|Gj-i|. Byinspection, for k 2 = k u (§|G fel _ 1 |+2 v /2|G fcl _ 1 |) 1 - 5 < 0.7|G fcl _x| 
provided |G fcl _i| > N Q = 10 5 . 

Assume the claim holds for k 2 — k\ = k — 1. 

fci+fc fei+i+fe-i 

X //, 1 < i^ fel r- 5 + E i^i 1 - 5 

j=fei j—^+i 

< |H fcl | 1 - B +0.7|G fcl | 1 - B 

< (6\G kl -i\ + 2v/2|G fcl _i|) 16 + 0.7 ((1 - 0)|G fel _x| + 2 v /2|G fcl _ 1 |) ^ 

where 9 is the balance parameter of the separator S kl ■ In the first inequality we 
have used the inductive assumption. Using the convexity of the above expres- 
sion, it can be bounded by setting 9 = g, which satisfies the lemma provided 
that all graphs have at least N = 10 5 nodes. QED 



Lemma 3.8 Assume T(n) < ain 15 logn for every No < n < \Gq\. Then, 
T(\H 1 \) + 2j:^ 2 T(\H ] \) <0.98a 1 |G | 1 - 5 log|G | 

Proof: Let 9\ denote the balance parameter for separator S\. 
T{\H l \) + 2Y,T{\H j \) < «i l^il 1 ' 5 log |G |+2ai log |G„| £ |ff/- 5 

< axl^il 1 ' 5 log |G | + 2 • 0.7a! log IGolldl 1 - 5 

< a 1 (0 1 |G o | + 2y2]G^) 1 - 5 log|G o | + 

2 • 0.7a! log |G | ((1 - 9i)\G \ + 2^2\G^j ^ , 

where in the second inequality we have used Lemma |3.7| Using the convexity 
of the above expression, it can be bounded by setting 9 = |, which yields the 
desired bound 0.98ai |G | 15 log |G | provided |G | > N . QED 



Lemma 3.9 for every i, every node of Ci belongs to some Jordan separator Sj. 
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Proof: Ci, the external face of Hi consists of nodes that either belong to Si 
or to a face of Gj_i that is not triangulated. To see that, consider a clockwise 
traversal of Si. For every two consecutive nodes of 5, that are connected in 
i by an edge, Si coincides with that edge (see Section 2.3), so it belongs to 



G 

Ci. The only parts of Ci and Si that do not coincide correspond to Si crossing 
some non-triangulated face / of Gj_i, say at nodes u and v. In those cases, 
d consists of the clockwise sub-path / between u and v. Since Go is triangu- 
lated in the first line of the algorithm and Line 13 every face of Gj_i that is 
not triangulated corresponds to regions that were strictly enclosed by previous 
separators, or more formally, to a union of the portions of Go that are strictly 
enclosed by some Jordan separators in {Sj : j < i}. Therefore, every node on 
these faces belongs to some Jordan separator Sj , which proves the lemma. QED 



We can now put together the pieces to prove the running time stated in 
Theorem [TU 



Proof: (Of running time in Theorem |1.1| ) We have already argued that 
the time required for all recursive calls is bounded by J2j ^Xl-^jl) + ^(I-^pC?) I)- 
The work done outside the recursive calls is dominated by the single-source 
single-sink flow computations in Line 16 Each of these computations takes 

uired for the non-recursive calls is 

this is 0(E,|^l|Go|log|G |). 



0(| Go | log | Go |) time. The overall time rec 
thus O 

Since |£j 
we have 



G 1 log | G 1 • By lemma 
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U 

is 0{y/\Gi-i |) and since the size of the G^s decreases exponentially, 



O ^|^||Go|log|G | < a 2 |Go|log|G | WlG,-! 



< c4|G r- 5 log|G | 

for some constants a 2, a' 2 - 

The overall running time is therefore bounded by 

T(\G \) < T(\Hj\) + T{\H p(]) \) + a 2 |G r- 5 log |G | 



Assume inductively T(n) < aijGol 1 ' 5 log |C 1 f° r some constant a%. Then, by 
Lemma |3~8l T(|G | < 0.98ai|G | 15 log |G | + a 2 |Go| 15 log |G | which at most 
011 1 Gq I log I Go I for appropriate choice of cti. QED 



4 Converting a maximum feasible preflow into 
a maximum flow 

In this section we describe a linear time algorithm that, given a feasible preflow 
in a planar graph, converts it into a feasible flow of the same value. This 
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algorithm can be used to convert the maximum preflow output by Algorithm [T] 
into a maximum flow. This section contains no novel ideas and is included for 
completeness. A similar procedure was used in [14j . but was not described in 
detail. 

First, use the technique of Kaplan and Nussbaum [TS] to make the preflow 
acyclic. The running time of this step is dominated by a shortest paths compu- 
tation in the dual of the residual graph. This can be done in O(nlogn) using 
Dijkstra, or in linear time using [llj . 

Let / denote the acyclic feasible maximum preflow in G. Let p(v) denote 
the net inflow of node v. Let D denote the DAG induced by arcs with f{d) > 0. 
Reverse every arc of D and compute a topological order on the nodes of D. The 
following algorithm pushes back flow from nodes with positive net inflow to the 
sources and runs in linear time. Upon termination, / is a feasible maximum 
flow. 



Algorithm 2 An algorithm that converts acyclic preflow p on a DAG D into a 
flow. 

1: for v G D in topological order do 

2: if v is not a sink then 

3: while p(v) > do 

4: let uv be a dart, where u comes after v in topological order and 

f(d) > 
5: x := mm{ f(d),p(v)} 

6: f(d) := f(d) - x 

7: p(v) :— p(v) — x 
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